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It has recently been argued by Copeland et al. that in eleven dimensions two 
orbifold planes can collide and bounce in a regular way, even when the bulk metric 
is perturbed away from Milne spacetime to a Kasner solution. In this paper, we 
f^ ' point out that as a consequence the global "phoenix" structure of the cyclic universe 

is significantly enriched. Spatially separated regions, with different density fluctu- 
ation amplitudes as well as different non-gaussian characteristics, are all physically 

CN ■ realized. Those regions containing by far the most structure are specified by a fluctu- 

CN ■ 

ation amplitude of Q ~ 10~'^'^ and local non-gaussianity parameters f^L ~ O(±10) 

and qml ~ C'(— 10^), in agreement with current observations. 
r^ ■ I. INTRODUCTION 
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In cosmology, it is tempting to declare that we should only ever think about what hap- 
pened in our past lightcone, as everything outside of it is inaccessible to observations. How- 
ever, several well-motivated models of the early universe predict the existence of regions 



f^^ . of the universe, with a whole range of physical properties, that are far separated in space 

and/or in time from our observable region of the universe. Hence it may be important to 
understand the characteristics of these far removed regions, at least on a theoretical level. 
It certainly seems of interest to explore all the consequences of a given cosmological model, 
whether these appear to be immediately observable or not; this is crucial in assessing the 
logic and internal consistency of the model in question, and hence crucial in determining 
whether a given cosmological model satisfies us intellectually. Such an understanding has 
concrete implications, as the perception of which properties of our region of the universe we 
consider to be determined by mathematical necessity rather than historical accident guides 
us in devising future research avenues. 

The question of additional universes is particularly acute for eternal inflation [1], since the 
number and variety of universes are both infinite^. Hence all possible universes are produced 



^ We will occasionally use the expressions "universe" and "region of the universe" interchangeably, when 
this does not seem confusing. 



an infinite number of times, so tliat, taken by itself, tlie model of eternal inflation is not 
predictive. In order to make it predictive one needs to add a measure, which regulates the 
infinities in question and leads to definite predictions. Unfortunately it turns out that all 
predictions depend sensitively on the measure/regulator used [2] - in other words, unlike 
the situation in QED for example, all predictions are regulator-dependent. Moreover, some 
of the simplest measures turn out to be in blatant conflict with observations. Hence we are 
certainly missing something important, and one may even conjecture that eternal inflation, 
as currently envisaged, is unphysical and does not occur. 

In this paper, we will examine the question of predictivity in an alternative cosmological 
model, namely the cyclic theory of the universe proposed by Steinhardt and Turok [3]. 
In doing so, we will extend existing studies of the global structure of the cyclic universe, 
previously described as the phoenix universe [4, 5]. Building on recent results by Copeland, 
Niz and Turok regarding brane collisions [6], we will show that the cyclic universe also leads 
to vastly separated regions with distinct physical properties. Hence one may be worried that 
a similar ambiguity as the one inherent to eternal inflation will come to plague the model. 
However, as we will describe, the cyclic model predicts strong correlations between a large 
number of physical observables in a given region, such as its spatial flatness, the amount 
of structure it contains and the non-gaussian statistics of its density perturbations. These 
correlations are entirely determined by the dynamics of the model, and thus they lead to clear, 
testable, predictions. For the speciflc model described below, it turns out that those regions 
that are flat and contain the most structure are characterized by a fluctuation amplitude of 
Q ~ lO^'^'^ and local non-gaussianity parameters fiy^ ~ O{±10) and g^L ~ 0(— 10"^). Thus 
the model is both in agreement with current observations, and falsifiable by upcoming CMB 
experiments. 

Before proceeding, we should clarify that the cyclic model contains a number of open 
questions. The most important such open question has to do with the bounce that connects 
the ekpyrotic contracting phase with the phase of expansion. This big crunch/big bang 
transition has been modeled by a collision of orbifold branes in M-theory [7]. To date, it 
has only been possible to analyze such a collision at the semi-classical level [6, 8], and not 
at the full quantum level. In the present work, we take the semi-classical results seriously, 
and explore their consequences. Thus, the present work is only applicable insofar as the 
semi-classical calculations accurately reflect the physics of the bounce. But, to the extent 
that our analysis yields promising results, it also motivates further work on the quantum 
treatment of the crunch/bang transition. 

The plan of the paper is as follows: we flrst provide a brief review of the cyclic universe. 



focussing on those aspects that are most relevant to the present paper. Section III then 
discusses the global structure of the cyclic universe, reviewing the picture as it was under- 
stood up to now, and extending it by using the results of recent semi-classical calculations 
of colliding orbifolds. The new global phoenix structure of the cyclic universe, together with 
its implications, is discussed. Conclusions and further directions are to be found in the 
discussions section. 

II. THE CYCLIC UNIVERSE 

In the cyclic universe [3, 9, 10], periods of expansion alternate with periods of contraction. 
During the expanding phase, the energy density of the universe is successively dominated 
by radiation, dark matter and dark energy. During this phase, the universe grows by a huge 
factor, of the order of exp{60 + A'^e), where A'^^^e denotes the number of e- folds of dark energy 
domination. Towards the end of the expanding phase, the scalar field responsible for dark 
energy starts to roll faster down its potential, until eventually the potential becomes negative 
- see Fig. 1. Shortly thereafter, the universe reverts from expansion to contraction, and an 
ekpyrotic phase ensues. The ekpyrotic phase is characterized by a very stiff equation of state 
{i.e. the pressure is much larger than the energy density during this phase), which has the 
consequence of suppressing both the average curvature and any anisotropies in the curvature 
of the universe, thus keeping/rendering the universe flat and isotropic [11]. During this 
phase, the scale factor of the universe shrinks by a very modest amount, while the Hubble 
rate increases tremendously. The ekpyrotic phase is followed by a brief phase of scalar 
kinetic energy domination, and during this phase the crunch/bang transition, which will be 
discussed in detail in the next section, occurs. Thus, the phases of expansion and contraction 
are highly asymmetric: in particular, over the course of one cycle, the universe grows by a 
huge amount. Only local quantities, such as the Hubble rate for example, behave cyclically, 
and return to the same value cycle after cycle. 

During the ekpyrotic phase, as the scalar field rolls down the steep potential, quantum 
fiuctuations are amplified and stretched beyond the horizon [12]. This is because, although 
the scale factor remains almost constant, the Hubble rate grows rapidly or, equivalently, the 
Hubble horizon shrinks rapidly. Thus fluctuation modes of successively smaller wavelengths 
find themselves outside of the horizon. In this way, the ekpyrotic phase generates classical 
perturbations, similarly to infiation. However, it turns out that the fiuctuations in the field 
driving the ekpyrotic phase have a blue spectrum, with spectral index n^ ^ 3, and a small 
amplitude [12-15]. Thus these perturbations are irrelevant on cosmological scales of interest. 
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FIG. 1: The cyclic universe potential and its main phases of evolution. 



To see how scale-invariant curvature perturbations are generated, it is instructive to first 
consider the embedding of the cychc model in string theory. The higher-dimensional view of 
the cyclic universe is set in heterotic M-theory, which describes 11-dimensional supergravity 
compactified on a 6-dimensional Calabi-Yau space and with one spatial dimension being a 
line segment (orbifold) S^ /'L2 [16, 17]. The (10-dimensional) endpoints of the line segment 
are called orbifold branes. For the quantum consistency of the theory {i.e. in order for 
anomalies to be absent), Hofava and Witten have shown that these branes need to have E^ 
gauge fields that reside on their worldvolume [18, 19]. These gauge groups are large enough 
to easily contain the Standard Model of particle physics, and explicit realizations of the (su- 
persymmetric) Standard Model have been realized in this context by adding supplementary 
branes in between the orbifold branes [20] (these supplementary branes do not concern us 
here). Thus, the 4-dimensional low energy effective theory of this setting can accommodate 
what we know about particle physics, while also containing additional scalar fields that can 
be relevant for cosmology. There are always at least two such scalars (moduli), one of which 
parameterizes the distance between the orbifold branes, and the other the volume of the 
Calabi-Yau space. Thus, from the higher- dimensional point of view, it is over-restrictive to 
consider a single scalar field. 

As we will review in section III, the 4-dimensional effective action for these 'universal' 
scalars is given by 
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where <^i and <^2 are related by a field redefinition to the inter-brane distance and the Calabi- 



Yau volume [21]. If we now assume that both scalars obtain an ekpyrotic-type potential^, 

V^(0i, 02) = -Vie-'-'^-' - V2e-'^'*>\ (2) 

where Vi, V2 are constants and ci, C2 slowly varying functions of 0i, 02, then both the back- 
ground dynamics and the characteristics of the generated perturbations change drastically. 
It turns out to be much more natural to discuss the dynamics in terms of the rotated scalars 
cr and s pointing transverse and perpendicular to the field velocity respectively [24, 25]; they 
are defined, up to unimportant additive constants which we will fix below, via 

0101 + 0202 0102 - 0201 .„^ 

a = : , S= : , (3) 

a a 

with a = (0^ + 02)^''^- It is also useful to define the angle 6 of the trajectory in field space, 

via [26] 

a ^^ ■ a ^"^ (A\ 

cosu = ^, smcy = ^. (4) 

a (7 

Then the potential can be re-expressed as 

V = -Voe^^-[l + es' + ^e'^'s' + ^e's' + ■■■], (5) 

where for exact exponentials of the form (2), one has k^ = 2v2(c^ — cl)/\ciC2\ and K4 = 
4(c5 + C2)/(c^C2(c^ + C2)). Here e, which is defined by 1/e = 2c^ + 2c2 and characterizes the 
steepness of the potential, is related to the equation of state w via e = |(l+w), and is assumed 
to be slowly varying. The potential describes an ekpyrotic direction along a, combined with 
a transverse unstable direction along s. This instability is key to what follows. It implies 
that there is a ridge at s = 0, and the evolution along this ridge is characterized by the 
ekpyrotic scaling solution 

a{t) = (-t)^/^ a = -J- In (-^/eVot) s = 0, (6) 

with the angle 9 being constant. 

As it is (7 that drives the ekpyrotic phase, this scalar field develops fiuctuations with 
a blue spectrum. However, as shown in [27], the transverse field s obtains (isocurvature) 
fiuctuations with an amplitude and spectral index given by 

Qs ~ |eK„,nr/' (7) 

n.-l = --^, (8) 



^ No such potential has been exphcitly derived from string theory to date, but there are indications that 
a potential of this form could arise due to membrane instantons stretching between the orbifold branes 

[22, 23]. 
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where the meaning of Kum will be explained momentarily. If e ~ (9(10^) or more, the 
spectrum is close to scale-invariant. In fact, given that the ekpyrotic phase must come to an 
end, the potential is expected to be less steep in the bottom half, leading to a natural range 
of 0.97 <n, < 1.02 [27]. 

But what we observe in the CMB are density fluctuations that are caused by fluctuations 
in the curvature of the universe. Thus, to match the model with data, we must calculate 
its predictions for the curvature perturbation (. On large scales, and at linear order, the 
curvature perturbation evolves according to [26] 

C = -—96s = J -96s. (9) 

Hence, if the trajectory in scalar fleld space turns {9 ^ 0), then the isocurvature pertur- 
bations 6s act as a source for the curvature perturbations (, which, given the absence of 
momentum-dependence in (9), will acquire the identical spectrum [28, 29]. Such a conver- 
sion of isocurvature into curvature perturbations can occur in two ways: 

• Ekpyrotic Conversion: if the trajectory is not localized close enough to the ridge during 
the ekpyrotic phase, then the trajectory will turn at a value of the potential V^um, 
and subsequently roll down one of the steep sides of the potential^ [24, 30, 31]. The 
amplitude of the nearly scale-invariant curvature perturbations that are generated in 
this process can be estimated by combining (7) and (9), yielding 

Ql ^ iKurnl- (10) 

Thus, the further down the potential the turn occurs, the more structure will be 
produced in the corresponding region of the universe. As shown in [32, 33], this 
conversion mechanism also generates large local non-gaussian corrections: writing the 
full curvature perturbation as a linear, gaussian piece (g plus correction terms as 

( = (g + lfNLCg + ^gNLC'g, (11) 

the predictions for the bi- and trispectrum are 

5 
/nl = ~T7:ci Ekpyrotic Conversion (12) 

9NL = —4, (13) 



^ During this further rolhng down, the universe is stih in the ekpyrotic phase, albeit with a shghtly different 
equation of state. However, now there is no further transverse instabihty, and the additional isocurvature 
perturbations that are produced during this stage of evolution have a blue spectrum and are irrelevant for 
cosmology, even if they should get converted into curvature perturbations at a later time due to a further 
turning of the trajectory. 
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where we have assumed that after the turn the trajectory evolves approximately in the 
02 direction (otherwise the formulae above should contain C2 instead of ci). 

• Kinetic Conversion: if the trajectory is localized close enough to the ridge during the 
ekpyrotic phase, so that the field makes it all the way down, then, as discussed in more 
detail in the next section, the trajectory automatically turns during the kinetic phase, 
when e ^ 3 [27]. In this case, the amplitude can be estimated again by combining (7) 
and (9) to be 

Qc~eek|Knin|, (14) 

where Knin stands for the potential value at the bottom of the ekpyrotic potential. 
Thus, compared to the maximal amplitude achieved by the ekpyrotic conversion mech- 
anism, the amplitude Qc_ is now enhanced by a factor of ^/t^ [34] . Detailed calculations 
show that for realistic values, e.g. Cck ~ 50, the enhancement in Q^ is only a factor 
of about 2, but nevertheless there is an enhancement. For this case of kinetic conver- 
sion, significant non-gaussian corrections are also produced, with detailed calculations 
yielding [35-38] 

3 

Jnl = 7T'^3\/e + 5 Kinetic Conversion (15) 

Qnl = (3/^4 + -^4 - 40) e, (16) 

where we may expect ks, /€4 ~ 0{1). 

Note that, as presented here and as conceived in the literature so far [24, 39], the two 
mechanisms of conversion appear to be rather distinct, and appear to depend on the specifics 
of a given model, in particular on the initial conditions assumed or generated at the beginning 
of the ekpyrotic phase. As we will demonstrate now, in the cyclic model of the universe this 
is not so and in fact all field trajectories are simultaneously physically realized! To make 
this statement explicit, we must look in more detail at the brane collision and the associated 
bounce phase. 

III. BRANE COLLISIONS AND THE PHOENIX UNIVERSE 

The central part of the cyclic universe is the bounce phase, where the transition from 
contraction to expansion occurs. In a fiat universe, such a transition can only occur if the 
null energy condition is violated, and thus we need physics that goes beyond the ordinary 
4-dimensional Standard Model of particle physics to describe such an event. Indeed, string 



theory contains many objects that violate the null energy condition (such as orientifolds, 
negative-tension orbifolds, anti-branes), but which nevertheless do not necessarily lead to 
inconsistencies regarding, say, the second law of thermodynamics. We will focus on the case 
where the bounce is caused by the collision of a positive- and a negative-tension orbifold 
plane in (heterotic) M-theory. Such a collision is singular from the 4-dimensional point of 
view, and thus the null energy condition is only violated at a single moment in time, and 
not over an extended period. There have also been attempts to model non-singular bounces 
[30, 31, 40], but it is not clear yet whether such descriptions are consistent or not [41, 42], and 
whether such non-singular bounces can arise in string theory or not [43-46] . The (classically 
singular) collision of orbifold planes currently represents the most promising description of 
a bounce. This brane collision has been analyzed both classically and semi-classically in a 
number of works [6, 8, 47-49], and we will base our study of the cyclic universe on these 
works. 

It was shown by Hofava and Witten that the strongly coupled heterotic string theory can 
be envisioned as 11- dimensional supergravity with one spatial direction being a line segment, 
or orbifold S'i/Z2 [18]. The size of this orbifold dimension determines the string coupling 
constant (which is larger, the larger the orbifold), and the endpoints of the line segment 
describe two orbifold planes of opposite tensions. For the quantum consistency of the theory, 
i. e. for the absence of anomalies, these two orbifold planes must contain gauge fields with 
symmetry group Eg. In order to obtain a realistic model of particle physics, 6 of the spatial 
dimensions should be compact and in the shape of a Calabi-Yau manifold (in order to ensure 
minimal supersymmetry in the 4-dimensional effective theory). It turns out that, in order 
to reproduce the observed value of Newton's constant, one needs to consider configurations 
in which the Calabi-Yau radius is smaller by one or two orders of magnitude than the 
orbifold dimension [50]. Thus it makes more sense to dimensionally reduce the original 
1 1-dimensional theory on a Calabi-Yau manifold first, and then investigate whether this 
theory has braneworld-type solutions, in which two branes are separated along an orbifold 
dimension. Such a dimensional reduction was performed in [16], with the metric being 
reduced according to 

dsl, = V~?'dsl + V^(,'dsly, (17) 

where Vcy denotes the volume of the Calabi-Yau manifold. The resulting 5-dimensional 



effective theory, called lieterotic M-tlieory, has action 
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R - \Vc^d^VcYd"'VcY - ^C?VcY 



^'^5 L J4d,H=+l i4d,«=-l J 



where we have kept only the cosmologically relevant fields, i.e. gravity and the Calabi-Yau 
volume modulus Vcy (the index m runs over time and 4 spatial dimensions). Here a denotes 
the amount of 4-form fiux wrapping a 4-cycle of the Calabi-Yau manifold, and must be non- 
zero in order to obtain braneworld solutions. In fact, as by assumption a 7^ 0, the vacuum 
of this theory is not given by 5-dimensional Minkowski space, but rather by the braneworld 
spacetime [16] 

ds" = h^/'>{y)[A^-dt' + dx^) + B^dy^], 
Vcy = Bh'"''{y), 
h{y) = 5ay + C, (19) 

where A, B and C are arbitrary constants. The y coordinate is taken to span the orbifold 
S^ /'L2 with fixed points at 1/ = ±1. In an 'extended' picture of the solution, obtained 
by Z2-refiecting the solution across the branes, there is a downward-pointing kink at y = 
— 1 (corresponding to a negative-tension brane) and an upward-pointing kink at y = +1 
(corresponding to a positive-tension brane). By inspection, one can see that the volume 
of the Calabi-Yau manifold and the distance between the boundary branes are determined 
in terms of the moduli B and C, while the scale factors on the branes are determined in 
terms of A and C. In a cosmological context, we are interested in describing the motion 
(and eventual collision) of these boundary branes. In order to do so, we can employ the 
moduli space approximation, which is valid as long as the branes move slowly (this will turn 
out to be a consistent requirement). The moduli space approximation consists of letting 
the 'constants' A, B, C become functions of time (for convenience, we use a conformal time 
coordinate r here). These moduli describe the lightest fields in the theory (as shown in 
[51], all other perturbations are massive) and the action specifying their dynamics can be 
obtained by plugging the now time-dependent moduli back in the action, yielding [21] 



4rf ^ vA/ 12\BJ AB 



1 U ^a 3 I-l^C 



-C'-'r-^ri^l (20) 
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where we have defined 

In = f dy h- = ,^^^A {C + 5a)("+i) - (C - ha)^-+% (21) 

This action can be greatly simphfied by introducing the field redefinitions 

a^ = A^Bh, (22) 

5 

g<^i/V2 _ B{Lif/\ (23) 

02 = -^ fdCih)-' [9(/„2)2 + 16/_iJ3]'/'. (24) 

^U y 5 >- 5 5 5-' 

The definition (24) can be inverted to give 

[(1 + g2v/2739i2N5/4 + (1 _ e^y273'^2)5/4l 

C = 5q; ^= ^= . (25) 

(1 + e^V2/3'^2\5/4 - (1 - e2V2/3'^2)5/4j 

In terms of a, 0i and 02? the moduli space action (20) then reduces to the remarkably simple 
form [21] 

^mod = / [-3a'' + la\<p[' + 0;^)]. (26) 

The minus sign in front of the kinetic term for a is characteristic of gravity, and in fact this 
is the action for gravity with scale factor a and two minimally coupled scalar fields that was 
used in the previous section, Eq. (1). The resulting equations of motion are solved by 

a" = a^(ro-r), (27) 

^cos^ln[0i,o(ro-r)], (28) 



-y|sin^ln[02,o(ro-r)], (29) 



where oq, 0, tq, 01^0, and 02,o are constants of integration. Thus, in the absence of a potential, 
the solutions to the equations of motion correspond to straight line trajectories at angle 6 
in the scalar field space spanned by 0i and 02. 

This scalar field space has one peculiar feature, which is best understood by first using 
the metrics (17) and (19) to relate 11-dimensional geometric quantities, in particular the 
inter-brane distance du and the Calabi-Yau volume Vcy± at the locations y = ±1, to the 
4-dimensional fields: 

VcY± = {2a)^^^e't'^/^ { ^cosh v/2/302 j ^ ^^^ ^ ^^^^^_^^^^^^ ^^^^ 

sinh a/2/302 
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FIG. 2: Examples of trajectories at different angles in scalar field space. The solid line trajectory 
corresponds to the kinetic conversion mechanism, and the dashed trajectories to the ekpyrotic 
conversion mechanism. 



The approximate relations on the right hand side exhibit the dependence of du and Vcy± 
on the scalar fields as the branes become close, i.e. as 02 — ^ — oo. Eq. (31) shows that the 
volume of the Calabi-Yau manifold at the location of the negative-tension boundary brane 
would become negative, if 02 became positive. As shown in [52], this unphysical region 
is generally not accessible, and the axis 02 = acts as a repulsive boundary - in other 
words, there is an effective repulsive potential near this axis. Thus trajectories approaching 
the boundary bend in its vicinity, and this turning automatically causes the conversion of 
isocurvature perturbations into curvature perturbations, as used above in the discussion of 
the kinetic conversion mechanism - see also Fig. 2 for an illustration. 

Fig. 2 depicts three different field trajectories, the solid line corresponding to a trajectory 
that makes it all the way down the ridge in the ekpyrotic potential and turns at the 02 = 
boundary before shooting off to — oo where the brane collision occurs. The dashed lines in 
the figure show two examples of trajectories that turn before reaching the bottom of the 
potential. These correspond to the ekpyrotic conversion mechanism. After their turn, they 
keep rolling down a steep side of the potential, and then shoot off to — oo as well, either 
directly, or after bouncing off the moduli space boundary at 02 = first. If such a second 
turn occurs, this will have little effect on the cosmological density perturbations, since the 
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additional isocurvature modes that are generated after the first turn have a blue spectrum, 
and have wavelengths that are too small to be of interest. In fact, the currently observable 
modes are those generated between about 50 and 60 e-folds before the end of the ekpyrotic 
phase, i.e. when the scalar field is a little more than halfway down the potential (assuming 
the potential minimum is approximately at the GUT scale). This means that for trajectories 
turning during the ekpyrotic phase, the turn must occur in the lower half of the potential if 
scale-invariant perturbations are to be generated on CMB scales. One may also wonder what 
happens to trajectories that turn early on. In order to conjecture about the fate of those 
trajectories, it is useful to re-consider the general requirements on the ekpyrotic potential: 
these are that it must be steep and negative, with a transverse unstable direction. It seems 
to be too much to assume that this transverse unstable direction will continue very far away 
from the ridge. It seems more reasonable to guess that trajectories that turn early will not 
undergo a long second ekpyrotic phase down the side of the potential, but rather that the 
ekpyrotic phase will come to a premature end. If this is so, the regions corresponding to 
those trajectories will develop strong curvature anisotropics in the approach to the brane 
collision, and will undergo a BKL-type mixmaster crunch [53]. We conjecture that these 
highly curved crunches will result in an over-production of particles at the brane collision, 
and that these regions will remain in a gravitationally collapsed state - locally, the branes 
will stick together. These regions stop growing and stop cycling. 

It is worth studying those trajectories that approach the brane collision after a long 
ekpyrotic phase in more detail. The collision occurs as (pi, 02 — t- — oo, i.e. the trajectories of 
interest are those described by Eqs. (27) to (29), with tt < 6 < ^. We can choose Tq = 0, 
so that the collision occurs as r — ?> 0^. Using (17) and (19), we can lift these solutions to 11 
dimensions, where (in the 02 — ^ ^oo limit) the metric becomes 

9^. ~ r/^,(2a)-i/^a2e-^'^^/(='^)+'^^/^ ~ (_^)i+5cose/(2V3)-sine/2 (33) 

gyy ^ (8«)-ie^'^i/(='^)+2-^^2/^ ~ (_^)-2cos9/^-2sine ^33) 

gab ~ (?,^,^(2a)l/4e'^^/(3v^)-'^^/^ ~ (_^)-cose/(2^)+sin9/2^ (34) 

with the a, b indices denoting the Calabi-Yau directions, and the intrinsic Calabi-Yau metric 
being denoted Qab'. Defining a new time variable T via 

_^ ^ rrp\4:/ (6+5 cose/ \/3-sme) /3g\ 

the metric can be brought into Kasner form 

ds^ = -dT^ + Y^ T^P^dxl (36) 

j=l,...,10 
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where i runs over all spatial dimensions. Then, denoting the exponents for the ordinary 
spatial dimensions by pi = p2 = ps = psd, the one for the orbifold by Py and those for the 
Calabi-Yau directions by P5 = Pe = Pr = Ps = Pd = Pw = Pcy, we get 

P^d = '^- —-Z T-75 — (37) 

+ 5 cos 6 /y/ 6 — sin 6 

-4cos^/y3-4sin^ 

6 + 5 cos 9/VS — sin 9 

- cos ^/a/3 + sin ^ .. 

^^^ = 7T~K — ^77f — ~^- ^^^^ 

+ 5 cos u/y/6 — sm u 

These solutions correspond to 11-dimensional Kasner solutions, with "^Pi = ^p'i = ^, but 

where the three ordinary spatial dimensions as well as the Calabi-Yau dimensions remain 

isotropic amongst themselves. In fact, as the orbifold dimension shrinks, the ordinary three 

spatial dimensions grow while the Calabi-Yau shrinks ii it < 9 < ^, and vice versa for ^ < 

9 < ^. The intermediate value 9 = ^ corresponds to the special case where the spacetime is 

Milne, with both the ordinary and the Calabi-Yau spatial dimensions approaching a constant 

value at the collision. 

The Milne solution was analyzed semi-classically in [8], and the Kasner solutions more 
recently in [6]. In all of these solutions, the curvature is small in the approach to the 
collision, and thus curvature corrections are suppressed. In fact, we may picture the effect 
of the preceding ekpyrotic phase as flattening and parallelizing the branes locally. As the 
branes become close, the lightest states in the theory correspond to M2-branes stretching 
between the two boundary branes. Their equations of motion remain regular at the collision, 
provided the Kasner solutions are "close" to Milne in the sense that Pi+Py/2 > for all i. This 
condition means that no direction should expand faster than (the square root of) the scale 
factor of the orbifold shrinks, and it is straightforward to verify that this condition is satisfied 
for all solutions above. Moreover, using semi-classical methods, the amount of particle {i.e. 
membrane) production at the collision was estimated and found to be proportional to the 
collision rapidity. For non-relativistic collision speeds, the amount of particle production 
is small, and we conjecture that universes undergoing such a bounce re-emerge unscathed 
and proceed to a new expanding phase. By contrast, for relativistic collision speeds, we 
may assume that there is so much particle production that the corresponding regions of 
the universe over-close and re-collapse rapidly. Such regions, along with those that did not 
undergo a prolonged ekpyrotic phase and enter the collision in a highly curved state, are 
filtered out by the brane collision, and effectively decouple from the expanding regions. 

We have just conjectured that there is a process of dynamical selection occurring at the 
big crunch/big bang transition. What are its consequences? Consider a region of space that 
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is at the end of its dark energy phase, and about to start ekpyrosis. As shown in [4], in 
scalar field space this region corresponds to a small blob sitting at the top of the ridge in the 
ekpyrotic potential, with a small spread As of field values along the transverse direction and 
a negligible spread Act along the ridge direction. Let A^tum denote the number of e-folds of 
ekpyrosis that occur before the trajectory turns. As discussed earlier, scale-invariant curva- 
ture perturbations on observable scales are only produced when this number is greater than 
about 70, and here we will assume this to be the case. Then a fraction e^^'""" of trajectories 
will roll down the ridge up to the potential height Vtum where the trajectory will turn, thus 
producing curvature perturbations with amplitude Q(^ ~ iKuml^^^ and non-gaussian correc- 
tions of the "ekpyrotic conversion" type. Afterwards, the trajectory will continue along one 
of the steep sides of the potential until eventually the potential bottoms out, becomes unim- 
portant, and a kinetic phase ensues. The trajectory will now shoot off to — oo in scalar field 
space, and there is no further instability along the rest of the trajectory. The brane collision 
will occur with a collision velocity yo determined by the minimum of the potential Vmm- As 
shown in [54], requiring this velocity to be non-relativistic implies a bound |Kiim| ^ 10^^*^ 
in Planck units. For now, let us assume that this bound is saturated. This implies that 
rolling down the entire ridge in the potential corresponds to 120 e-folds of ekpyrosis. Thus, 
a fraction e^^^" of all trajectories undergo the full ekpyrotic phase and turn only during the 
kinetic phase. These "kinetic conversion" trajectories generate curvature perturbations with 
the largest amplitude, Q^^ ~ iKiiinl^^^ ~ 10~'^'^ [54], with the corresponding non-gaussian 
corrections, and they lead to a brane collision with approximately the same velocity yo. Af- 
ter the brane collision, all regions that made it through will expand by the huge factor of 
g60+7Vde^ as described in section II. 

Above, we argued that it is unreasonable to assume that the turn in the trajectory can 
occur arbitrarily early during the ekpyrotic phase and still make it through the bounce. If 
the earliest that the turn can occur corresponds to Aturn,min, then we must compare the losses 
occurring because of the instability of the ekpyrotic phase with the amplification occurring 
in the expanding phase, i.e. we need Aturn,min < Nde + 60, or Nde > Atum.min — 60, in order 
for the cycling to be sustainable. Thus dark energy remains crucial in ensuring the survival 
of the cyclic universe, but since Aturn,min < 120, the constraint on the total amount of dark 
energy expansion is weakened from the earlier bound, which was N^e > 60. 

The resulting global structure is an extension of the phoenix universe of [4], see Fig. 3: 
during the ekpyrotic phase, most field trajectories turn too early and end their ekpyrotic 
fiattening phase prematurely. These regions end up in gravitationally collapsed structures, 
e.g. black holes. But some field trajectories, approximately e^^'""^"™'" in number, undergo 
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FIG. 3: A spatial slice through the phoenix universe. 



a long ekpyrotic phase, generate nearly scale-invariant curvature perturbations and make it 
through the bounce. These flat and potentially habitable regions get amplified hugely during 
the phases of radiation, matter and dark energy expansion, and they eventually spawn new 
cycles of evolution. In previous work [4], there were two assumptions that have now been 
relaxed: the first was that only the Milne bounce, i.e. only the special trajectory at angle 



e 



In 



, was viable. And the second was that this trajectory was assumed to correspond to 



the "kinetic conversion" mechanism. Now we see that both assumptions are unnecessary: 
the orientation of the ridge can be at any angle in scalar field space, and all trajectories that 
experience sufficient ekpyrosis make it through the bounce, as long as their angle is in the 
range n < 9 < ^. The consequence is that we now get a diverse structure in the phoenix 
universe, with a range of physical properties. The main segregation is still between large, 
fiat regions and tiny, highly curved ones. But within the large, fiat regions there is now 
considerable variety. 

Of the fiat regions, the rarest are those that make it all the way down the ekpyrotic 
potential and convert their isocurvature into curvature perturbations during the kinetic 
phase - but these regions also have the largest perturbation amplitude. By contrast, the 
regions that correspond to earlier turns of the trajectory will be far more numerous, but will 
also contain a lot less structure. Hence it is of interest to estimate the amount of structure 
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formed in these various large regions. The full ekpyrotic phase corresponds to about 

iVek = ^ In %^ ^ 120 (40) 

e-folds, where Vq denotes the present energy density of dark energy. The regions that turn 
at Kuril undergo only 

iVturn = ^ln^ (41) 

e- folds of ekpyrosis. These regions are thus more numerous by a factor e^-^k^^tum^ ^^^ corre- 
spondingly, they occupy a volume that is larger by a factor of 

The maximal density perturbation amplitude is achieved by the kinetic conversion regions - 
we will denote this amplitude Qmax- As discussed at the end of section II, for realistic values 
of the equation of state, this amplitude is about twice as large as the largest amplitude 
that can be achieved with the ekpyrotic conversion mechanism. Given that Q oc |l^n^, 
this implies that for a trajectory that turns at Kurn, the associated density perturbation 
amplitude is 

Qturn ^ ^iY^fl\ (43) 

If we assume the Press-Schechter model of structure formation [55], then in a given region of 
space, the probability for having a fluctuation with amplitude b is proportional to A^^^ ^^ ■ 
Hence, the number of fluctuations of amphtude b in the regions turning at Kurn, compared 
to the kinetic conversion regions, is 

''\ (44) 

The exponential suppression dominates over the power-law increase in volume, showing that 
there are far fewer large fluctuations in the "ekpyrotic conversion" regions. In fact, one 
can integrate over Vturn to estimate the number of fluctuations of amplitude b ~ Qmax in 
all ekpyrotic conversion regions compared to the kinetic conversion ones, obtaining (with 

(XX ^ upturn/ •'minj 

C^^ 2(^)2e-(^^'^)^ = /' \e'-'l-dx ^ 0.025. (45) 

Hence, there are about 40 times as many density fluctuations of magnitude 10~^ in the 
kinetic conversion regions than in all other regions combined! 



^ Vturn ^ Qturn 

1 ^-A2/Q2^^_ 


= 2(^) V 


^ /"/[ ^min 


Qmax ^ ' "^'turn 


Wmax 
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IV. DISCUSSION 

The best understood version of the cychc universe involves an ekpyrotic phase during 
which the potential is unstable. Because of this instability, nearly scale-invariant isocurva- 
ture perturbations are generated, which get converted into curvature perturbations as the 
trajectory in field space turns. This turn can happen in essentially two ways, either during 
the ekpyrotic phase or during the subsequent kinetic phase. It was thought up to now that 
these possibilities correspond to different cosmological models, depending on how the initial 
conditions are arranged. Here, we have shown that in the cyclic universe, where each cycle 
generates the "initial conditions" for the next one, all of these possibilities are physically 
realized in parallel. 

This leads to a global structure of the phoenix type, where large, flat regions now come 
in a variety of physical characteristics, but, interestingly, with strong correlations between 
a number of cosmological observables. The regions containing by far the most structure 
have a primordial density fluctuation amplitude of Qmax ~ 10~^'^, with local non-gaussian 
corrections specified by j'nl ~ O(±10) and qnl ~ (9(— 1000). There are far more numerous 
regions with a continuous range of ever smaller density perturbation amplitudes (starting at 
about Qmax/2) and local non-gaussianity parameters /nl ~ O{—10) and qnl ~ O(+1000). 
Even when counted together, these regions contain far fewer large amplitude fluctuations 
though than the regions having the maximal possible density fluctuation amplitude, and 
thus most galaxies are to be found in the Qmax regions. Also, as is typical in cyclic models, 
none of these regions are expected to contain large amplitude primordial gravitational waves. 
Finally, there are regions that are in a gravitationally collapsed state. Those are regions 
where the ekpyrotic phase ended prematurely, and which (we conjecture) underwent a high- 
curvature brane collision from which they emerged with an overproduction of matter, causing 
them to re-collapse rapidly. 

One can summarize the essential steps that lead to the global phoenix structure as follows: 
diversify, filter, amplify. The unstable potential diversifies the possible physical properties 
that the model produces, the brane collision filters out those regions that are sufficiently 
fiat, and the radiation, matter and dark energy phases amplify those regions that made 
it through the bounce in good shape. This sequence of steps contributes to making the 
model highly predictive. Another aspect enhancing its predictivity is the fact that during 
the longest smoothing phase, i.e. during the dark energy phase, the Hubble rate is lower 
than during the principal phases of structure formation {i.e. the ekpyrotic, radiation and 
matter phases). This implies that even large quantum fiuctuations during the dark energy 
phase do not change the structure of the model in an essential way (in contrast to eternal 



inflation). The only consequence of such fluctuations are that the corresponding regions 
expand a little more, and enter their next cycle a little later than the surrounding regions. 
A more quantitative analysis of this particular point will be left for future work. 

From a model-building perspective, the new phoenix universe presented here has two 
advantages over the "old" version. First, it shows that no fine-tuning is needed regarding 
the orientation of the unstable ridge of the ekpyrotic potential. And secondly, the constraint 
on the number of e-folds of dark energy needed to sustain the cyclic universe is weakened. 
In fact, in the new picture, it is likely that our predecessor-universe was a region with low 
Q and very little structure. 

In this framework, there are many opportunities for further study: one may wonder for 
example how sharp the boundaries between the different regions in the phoenix universe are, 
and whether these boundaries lead to any interesting effects. Evidently, it would be desirable 
to have a fully quantum calculation of the bounce, but such a treatment may still be quite 
far into the future. In the meantime, it would certainly be worthwhile to study more detailed 
models of the bounce, and perhaps investigate under what circumstances physical quantities 
can change. An interesting possibility is that the bounce may be a strong filter for many 
other physical quantities too, and thus may select a whole range of features of our universe 
that otherwise might appear to be merely historical accidents. 
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